Quantum crystal phase of electron liquid 
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A new bosonic, excitonic method for interacting electrons is developed. For two-dimensional 
electron liquid, it reveals a noncompensated quantum crystal phase ranging from the high density, 
or vanishing Coulomb interaction, limit: r s — 0, to a lower density r s ~ 10, where a condensation 
of exciton pairs and the loss of the Fermi step is possible. The correlation energy and the static 
susceptibility, calculated using the method, are in agreement with the independent studies. A proof 
of the instability of the Fermi surface in the limit r 3 — > is given. The concept of quantum crystal 
is capable to explain the metal-insulator transition in Si field-effect heterostructures. 
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Electron liquid is a system of non-relativistic electrons 
with Coulomb interaction Its quantum ground state 
depends on the density of electrons n only, or on the nor- 
malized inter-electron distance r s , with r s = e 2 m/ ' \Jisn in 
two dimensions (2D). Landau theory of Fermi liquid [3,0| 
and quantum Monte Carlo (QMC) method 0, 0| predict 
that paramagnetic, uniform liquid, with the same type of 
excitations as in the free Fermi gas, is stable in 2D from 
r s = to up to r s ~ 30, whereas the experiment in 2D 
Si field-effect heterostructures finds a metal-insulator 
transition (MIT) at r s ~ 9. Landau theory describes the 
Fermi liquid in terms of the two-particle scattering vertex 
r Q. The translational, Galilean and gauge symmetries 
of the Fermi liquid impose constraints on the vertex T 
Q, however, they alone are insufficient to determine it. 
While r is unknown for arbitrary r s the conditions for 
the stability of the Fermi liquid is expressed in terms of 
T 0. The QMC method applied to Fermi liquid or 
Wigner crystal is limited by the use of variational wave 
function and, therefore, studies the phase predetermined 
by the wave function guess. The QMC method was in- 
capable to reveal the instabilities of the Fermi liquid. 

In this paper I propose a new method for uniform elec- 
tron liquid (eliquid) where degrees of freedom are bosons 
representing excitons, the electron-hole pairs. The calcu- 
lated correlation energy of 2D electron liquid as a func- 
tion of r s agrees with the QMC results The method 
reproduces the Coulomb screening of the static potential 
[H|. New features found are two instabilities. First, the 
susceptibility on the Kohn momentum, twice the Fermi 
momentum pp, diverges at r s ~ 2 and this instability is 
independently proven to develop for any r s > 0. Second, 
the method reveals a Bose condensation of exciton pairs 
around r s ~ 10, accompanied by the loss of the Fermi 
step discontinuity. Qualitative interpretation of these 
findings is that charge density waves develops opening 
gaps on the Fermi circle. The ground state of the 2D 
eliquid is the quantum crystal [8| with, perhaps, trian- 
gular lattice and with the number of electrons that does 
not match the number of the lattice cites. Noncompen- 
sated electrons form two smaller, closed, hole-like Fermi 



circles. MIT is a point where the Fermi circles disappear 
and the quantum crystal transforms into the compen- 
sated Wigner crystal. 

Only 2D spin polarized eliquid is considered here. The 
ground state of the corresponding Fermi gas is: \G) = 
Yl ip + (s) |0). Everywhere, s and p stands for a momen- 
tum inside and outside the Fermi circle. Accordingly, the 
Hamiltonian is the sum of i) the Hartree-Fock (HF) term: 

Ho = E HF + Y £ hf(p)^^p ~ Y e H F{s)^t (!) 



where Ehf is the total HF energy, €hf(p) is the HF 
dispersion [![, and the Coulomb terms: ii) creating or 
destroying two excitons: 

Si<S 2 ,pi<P2 

(2) 

where J7(pip 2 ;s 2 si) = (E7(pi - si) - U(pi - s 2 )) <5(pi + 
P2 Si— Sa), t/(q) = 27re 2 /|q|„4 is the Coulomb potential, 
A is the area filled by eliquid, iii) scattering quasiparticles 
while preserving the number of excitons: 

h 22 = Y c/(piP2;p' 2 p'i)^ 1 ^ 2 Vv 2 V 1 + 

Pl<P2,P'i<P' 2 

+ Y C / (sis 2 ;s' 2 s'i)V's 1 V's 2 V'X^i ~ 

si<s 2 ,s'i<s' 2 

- Y c/ (pi s 2; s iP2)V'p 1 V' S iV' s + 2 '/'p 2 (3) 

pip 2 sis 2 

and iv) creating or destroying one exciton: 

H 12 = Y U(PS 3 ;S 2 S 1 ) (V^aa^Vii + h.C.) + 
ps 3 ,si<s 2 

Y E/(p 3 s;P2Pi) (ip+ 3 ip+ip P2 ip Pl + h.c.) (4) 



SP3,P1<P 2 



All quantum states that emerge during the temporal evo- 
lution of eliquid with the Hamiltonian {T][4|) are neutral 
and translationally invariant, i.e. with equal number of 
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quasiparticles and quasiholes and with zero total exci- 
tation momentum. Among these excited states there is 
a set of irreducible states with the property that quasi- 
particles could not be partitioned into two neutral sub- 
sets with both having the zero total momentum. Let 
(n, pi...p„, Si...s„) or (n,i), with zero total momentum, 
be an index of the irreducible n-exciton state and 

n n 

^ni = n<-IK ( 5 ) 

be the operator creating this state. They all commute: 
lKvA+ j }- = (6) 

A state of eliquid is partitioned into the irreducible states 
very much like a number is a product of the prime num- 
bers. In general this partition is not unique. For example, 
a four-exciton state with two excitons having momentum 
q and the other two excitons having momentum — q can 
be reduced into the irreducible two-exciton states in two 
different ways as shown in Fig.lA. 

Consider a bosonic state specified by the index (n, i) 
and the corresponding Bose operator a+ . We associate 
with the action of operators A in the electron Fock space 
the action of operators A in the boson Fock space: 

\exc)=A+ \G) <- |0..1 ni ..0> = A+ |0) (7) 

The boson Fock space is larger than the electron Fock 
space. First, the occupation number of any boson state 
can not exceed one. Second, a double excited state 
|0..1„i..l mj -..0) = A^A^j \G), which is correctly defined 
because of Eq.©, vanishes if at least one momentum en- 
ters both sets (ni) and (raj). Let f2 p be a set of the 
irreducible states (n, i) that include the momentum p. 
Then, the operator: 

n 

Ki ^PL.p^.^n^A (8) 

conforms to the electron Fock space limits if the number 
of slave bosons b p s equals to one in the ground state |G): 

b Pp+ E a ™«™ = 1 ( 9 ) 
(n,i)en p 

This representation is similar to the Schwinger boson 
representation for a spin. The analog of the Holstein- 
Primakoff representation of a spin in our case reads: 

j v mien Pj y m/eo Sj . 

(10) 

To avoid the double counting of the states in Fig.lA we 
exclude them in the a-boson space while a new boson: 
c 4m' wnere (4m) is (2«i) (2z 2 ) or (2ji)(2j 2 ) identical be- 
cause the momenta in these two sets are the same, will 
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FIG. 1: A. Insert. The four-exciton state that can be reduced 
in the two different ways: (2ii), (212) or (2,72). 
B. Static susceptibility in units of m/2ir vs r s , for N = 37 
and two wavevectors: the minimum on the lattice Q ~ 0.3pf 
and the nearest to the Kohn momentum: Q ~ 1.86pf- 

represent such a state. Let f2 q be a set of two-exciton 
states (2i), with the momentum q in one of the two 
particle- hole channels || or x, and let / q be the corre- 
sponding slave boson. Then, the operator 

A-ti = a piP2,s 1 s 2 ^Pl^Sl^P2^s 2 /q|| /q x (H) 

and the constraint: /+ / q + X)(2i)eo q a ti a ^i — 1' ex " 
elude the double-counting states. This constraint can 
be resolved to give the Holstein-Primakoff representa- 
tion. Namely, it is Eq. (flTj)l multiplied by two squareroots: 

0||,x yl -E(2m)en q a 2m a 2m, corresponding to the two 
double-counting channels. Representation of the higher 
many-exciton states may be more involved then the two- 
exciton state Eq. fTTT) . 

Our aim is to map the Hamiltonian (jTJUJ) into the boson 
space. We do this by expanding the Hamiltonian order 
by order in number of excitons. This expansion is similar 
to 1/5* expansion in the theory of antiferromagnet. It has 
been noted that 1/S expansion could be convergent 0], 
with the higher order terms contributing lesser correc- 
tions to the lower order terms, and this is due, perhaps, 
to the compactness of the spin manifold. Usual in the 
theory of Fermi liquid is the expansion in powers of r s 
[l|, and this, perhaps, is divergent, with long series and 
exact pattern of divergence required in order to make a 
useful resummation, like Borel in the theory of critical 
indices of the second order phase transition. 

The boson Hamiltonian is expanded in terms ordered 
by the exciton number. First is the two-exciton term: 

H 2= E ^(Pi,P2;s 2 Si) (a+ P2S2Sl +a PlP2S2Sl ) 

pi<P2,Sl<S 2 

(12) 
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TABLE I: Correlation energy per electron in atomic units and 
the Fermi step discontinuity 



r s 


1 


2 


3 


4 


5 




0.0235 


0.0194 


0.0167 


0.0149 


0.0135 


An{p F ) 


0.95 


0.88 


0.80 


0.72 


0.62 



and next is the four-exciton term 
H 4 = J2 a 

ij 



(13) 



Scattering takes place in 6 different channels: particle- 
particle (first line of Eq.Q), hole- hole (second line of 
Eq.Q), and four different particle-hole channels (third 
line of Eq.©). H 22 matrix includes the fermionic sign 
of the momenta ordering before and after scattering. 

Minimum of the Hamiltonian H2 + Hi with respect 
to a gives the correlation energy of eliquid in the lowest 
order. Quantization of momentum on the triangular lat- 
tice is used. For number of electrons N = 31, 37, 55, 61 
the ground state possesses D$h symmetry which reduces 
the dimension of the Fock space by twelve. The mo- 
mentum p is restricted by a maximum value Rpf, with 
typically R = 2.. 5. The dependence of E c on R is smooth 
and fitted by the law: E C (N,R) = E C (N) - k/R 13 , with 
(3 = 5. .8 for r s = 0..5. E C (N) as a function of N shows 
wide 'mesoscopic' variations, related to the electrostatic 
energy of N electrons put in the box with the size de- 
termined by the lattice spacing. E C (N) variations are 
correlated with the variations of the particle-hole gap on 
the Fermi circle A(iV), with A (00) = 0. In the limit 
r s «lwe find E C (N) = E c0 (N) - aA(N), for the quasi- 
particle dispersion: e = yj '£ 2 + A 2 . Therefore, A(iV) 
is also calculated and the best a is found to lessen the 
E C (N) variations. Then, E c q(N) is fitted by the law re- 
flecting the Coulomb singularity: E c o(N) = E c + ~//y/N. 
E c is shown in Table I and it is within 10% of the 
QMC result 0]. The next approximation to the bosonic 
Hamiltonian includes: H$ term that scatters two excitons 
into three-exciton states and Hq term that scatters be- 
tween three-exciton states as well as accounts for nonlin- 
ear three two-excitons processes. The occupation num- 
bers (if>pif)p) = YlieQ ( a 2i°2i)/2 and the corresponding 
Fermi step discontinuity have been found, see Table I. 
The Fermi step discontinuity vanishes at r s close to the 
point of the two-exciton Bose condensation. 

The lowest eigenvalue of the matrix H 22 is shown in 
Fig. 2. A bound state of two excitons with negative en- 
ergy occurs at r s > 9. .13 with 'mesoscopic' variations 
depending on N and excitons condense into this bound 
state. The configuration of two excitons, that contributes 
the most to the condensate wave function, has two holes 
with si = — S2 and two particles with pi = — P2. 

The response of eliquid to an external static poten- 
tial 2Vq cos(qr) can be also found using the boson repre- 




FIG. 2: The lowest energy in atomic units of two excitons as a 
function of r s for TV = 61, 55, 85, 37 in the order of increasing 
r s of the condensation point 



sentation. The potential breaks down the translational 
symmetry and, accordingly, an additional q-boson space 
with one- and two-exciton operators auq and d2i q is in- 
troduced. For small Vq, q-states are slightly excited and 
the nonlinearity is negligible. To exclude the double 

with two excitons of 
,4- 



counting of the states a^a^ lO), 

the state 2j having momenta ±q, the operators an 



modified as: 



'2j 



1 



The q-Hamiltonian 



for the eliquid response in the four-exciton order consists 
of the Coulomb and the potential terms: 



E 



a liq^jjq a ljq 



21 

ijq U ljq 



h.c 



Vr 



E 



+ rr22 



«- q + Oliq) 



+^qE( a 2j a Hq 



"j.q'V;C|"-', 



h.C.)yU) 



where <5 q is the momentum shift by q. Integrating out the 
variables an from the Hamiltonian H2 + H4 + H* we find 
an effective q-Hamiltonian in the q-space. Its minimum 
with respect to auq and a2i q gives the energy gain: Eq = 
— x(q)V 2 . The susceptibility x(l) as a function of r s 
is shown in Fig. IB. In the long range limit q <g. pf it, 
the lower curve, is suppressed by the Coulomb screening: 
x{l) ~ l/ r sPF, whereas the Kohn susceptibility x(%Pf), 
the upper curve, diverges at r s ~ 2. This divergence is 
consistent with the unrestricted HF study of 2D eliquid 
where larger N has been accessed and the instability 
towards the Wigner crystal has been found at r s ~ 1.4. 

We prove that the Fermi gas ground state \G), assumed 
in the perturbative studies of the eliquid [l|, 111] , is un- 
stable even in the limit r s < 1. Consider a general HF 
state: \Gjj) = Il s J2 P V'p Up S \0), created by the unitary 
rotation U in the momentum space. It has the energy: 
E HF = (l/2m)Tr(p 2 p) + 

+ \ E U (<ti t Tr (M Tr ( S -iP) - Tr (V<UP)] (15) 
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where p = UNU + is the density matrix, N is the di- 
agonal occupation matrix: N s — 1 and N p = 0. The 
parameterization: 

where z is the rotation matrix from s to p space, allows 
one to expand the energy Eq. (fT"5f in the vicinity of \G): 

Ehf = ^(e ff F(p) - £ffF(s))z s + p Zp S + 

ps 

+ X)(C(q) - ^(S - sOXs+q^s'+q.s' + 
qss' 

J2 (U(q) - U(s -s' + q))« s+q 4,s'- q + h.c) (17) 

q>0,ss' 

The instability momentum is |Q| = 2pi? and the corre- 
sponding exciton wave function is tjj(s) = z a +Q, 8 VA. For 
small deviations from the Fermi circle: r = pf — \s\ <C pf, 
the wave function factorizes: ip(s) — ip(r)ip(9), where 6 
is the angle along the Fermi circle. The function ij)(9) is 
positive within limits 6 max — Rq/pf where Rq = ^pftq 
is the largest and ro is the smallest momentum scales of 
the function tp(r). Then, the energy Eq. (fP7j) becomes: 

^ = S/T r\^- 9 -£drdr'r(rmr') log 

(18) 

where g = e 2 /(nvp) < 1 is the coupling. With the 
logarithmic accuracy the wave function depends on £ = 
log(r/ro). Provided yfgL ~ 1, where L = log(ilo/Vo), the 
HF term in the kinetic energy is small ~ gL. Variation of 
the functional Eq. fTg)) gives the equation for ijj(r) with 
the solution: tp(r) — acos( y /gS : )/r, for the eigenvalue 
^fgh — 7r/2 and within the range ro < r < i?o, where 
a ~ y/ro is the normalization. The energy Eq. (|18|) for 
this solution: Ehf = 2vfo? sin (2y/gL) /y/g, becomes 
negative at 2^/gL = ir and is small Ehf ~ —VfTq, for 
r s « 1, with the instability momentum scale: 

r w p F exp {~TT^/2Ti/r^j (19) 

This proof can be extended to 3D and to eliquids with the 
spin. The concept of the 'normal' Fermi liquid does not 
apply to electron liquids, it is rather a quantum crystal. 

In conclusion a new method for interacting electrons 
has been developed. It reveals that 2D electron liquid 
is a noncompensated quantum crystal at r s < 10 and 
is a fluid. Even the Wigner crystal, in the presence of 
pinning centers with the density rii, flows. Indeed, a uni- 
form force deforms the lattice and the stress builds up 
around the pinning cite making the quantum tunneling 
action S(r s ) to be different for the lattice slip along or 
counter the force. The conductivity of the Wigner crys- 
tal is <ic ~ (e 2 / Ti) yj n / rii exp(—S(r s )). The excitations of 



the eliquid falls into two types: the usual Landau quasi- 
particles in the vicinity of the Fermi surfaces and new 
Bose modes induced by the crystal order. The conduc- 
tivity of the noncompensated quantum crystal includes 
the quasiparticle and the crystal terms: a = op + oc- 
The classical Hall effect of the noncompensated quan- 
tum crystal is <jh = enc/H. Indeed, only the Lorentz 
and impurity forces are applied to electrons, but the lat- 
ter is a viscous force with no component perpendicular 
to the current. The Hall effect refines into the negative 
part of the hole- like Fermi circles (involving the umklapp 
processes) and the crystal part from the whole Brillouin 
zone given by the Chern number. MIT in this picture 
is a disappearance of the Fermi circles as the crystal or- 
der parameter and the particle-hole gaps grow. In the 
MIT point the quantum crystal transforms into a regu- 
lar Wigner crystal. For large crystal action S(r s ) ~ y/r^ 
the conductivity, a w ap, will show the activated be- 
haviour on the insulator side whereas on the metal side, 
where the charged carrier density n(T) is low, it will be 
governed by the 2D classical plasma law at finite temper- 
ature: <jl ~ (e 2 /ft)(e 2 -\/n(T)/T). In the parallel mag- 
netic field electrons in the crystal phase are easily polar- 
ized and all the tunneling probabilities, determining the 
quasiparticle mass and mobility, will decrease due to the 
Fermi sign repulsion. 
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